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We present a fluctuation relation for heat dissipation in a nonequilibrium system. A nonequilib- 
rium work is known to obey the fluctuation theorem in any time interval t. A heat, which differs 
from a work by an energy change, is shown to satisfy a modified fluctuation relation. Modification 
is brought by correlation between a heat and an energy change during nonequilibrium processes 
whose effect may not be negligible even in the t — > oo limit. The fluctuation relation is derived for 
overdamped Langevin equation systems, and tested in a linear diffusion system. 
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Fluctuations of thermodynamic quantities of nonequi- 
librium systems obey a universal relation referred to as 
fluctuation theorem (FT) [J-Q. Discovery of the FT 
leads to a great advance in nonequilibrium statistical me- 
chanics. Based on the FT, one can generalize the fluctu- 
ation dissipation relation to nonequilibrium systems fiol — 
12j and figure out fluctuations observed in experimental 
small-sized systems 13 - TB| . 



The FT for a quantity 1Z over a time interval t takes 
the form (e _TC ) = 1, where the average (•) is taken over 
a probability distribution for an initial state and over 
all time trajectories. Some quantities further satisfy the 
FT in the form P r {R)/P r {-R) = e R where P r {R) = 
(5(71 — R)) is a probability density function (PDF) for 
a nonequilibrium process and P r (R) for a corresponding 
reverse process. The latter is called the detailed FT and 
implies the former called the integral FT. 

Consider a system being in thermal equilibrium with 
a heat reservoir. We will set the temperature and the 
Boltzmann constant to unity. The system is driven into a 
nonequilibrium state if one adds a nonconservative force 
or applies a time-dependent perturbation. Then, there 
exist nonzero net flows of a nonequilibrium work W into 
the system and a heat dissipation Q into the reservoir. It 
is well established that the work W over a time interval t 
obeys the FT 0, Q . In addition, the total entropy change 
AStot = AS sys + A<S res with the system (reservoir) en- 
tropy S sys (S res ) satisfies the integral FT for an arbitrary 
initial state, and even the detailed FT for a steady state 
initial condition Q . Thermodynamic quantities are mea- 
surable experimentally from time trajectories in classical 
systems [l6[, while their experimental measurability in 
quantum systems is still an open issue [17 1. 

Fluctuations of heat Q, or entropy production AS res = 
Q/T, has also been attracting much interest [l8l - l26| . 
Note that a heat differs from a work by an energy change 
AS as Q = W — AS. When t becomes large, the system 
will reach a steady state with constant work and heat 
production rates on average. Hence one may expect the 
FT for heat in the large t limit where an energy change 
can be negligible (Q ~ W > AS). In fact, the FT for 



the heat production rate (Q/i) is derived formally in the 
t — > oo limit 0, Q • On the other hand, some model stud- 
ies demonstrate the FT for heat [25[ or failure of the FT 
in the t — > oo limit [ljl 0, 24 1. So, it is interesting to 



understand how and why the FT is violated for finite t 
and whether it is restored in the large t limit 21 , 22j ■ 



In this Letter, we present a fluctuation relation for 
heat, given in Eq. For any process, fluctuations are 
constrained by the energy conservation Q = W — AS. So 
a correlation between thermodynamic quantities plays an 
important role in characterizing the heat fluctuation. We 
find that the heat distribution satisfies a modified fluc- 
tuation relation that differs from the ordinary FT by a 
factor reflecting such a correlation. The fluctuation rela- 
tion is confirmed for a linear diffusion system analytically 
and numerically. Our work provides an insight into ori- 
gin for failure of the FT for heat for finite-i interval and 
possibly for infinite-i interval. 

We consider a dynamical system described by an over- 
damped Langevin equation 



dq(t) 
dt 



(1) 



where q = (qi,q2,'" ,qd) T is a configuration vector, 
/(<?) = (h(q)j2(q),--- ,fd(q)) T is a force, and £ = 
(Cij^2, • • • j£d) T is a white noise with 



(&(*)> =o, &(*)&(f)> = 2M(*-0 



(2) 



A damping coefficient and a noise strength are set to 
unity by rescaling t and q properly The force can be de- 
composed as / = f c + f nc , where f c (q) = - V q $(q) is a 
conservative force with a scalar potential energy function 
<f>(q) and f nc (q) is a nonconservative force. In this work, 
we focus on systems with a time-independent potential. 
We assume that the system is in thermal equilibrium fol- 
lowing the Boltzmann distribution P eq (q) oc e - *^ 9 ^ ini- 
tially at t = 0. Then it evolves into a nonequilibrium 
state due to the nonconservative force. 

When the system follows a path q(r) for a time interval 
< t < t, a nonequilibrium work done by the noncon- 
servative force, a heat dissipation, and an energy change 
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are given by functional W[q(r)] = J Q drq(r) • /„ c (<?(t)), 
Q[q(r)} = f* drq(r) ■ f(a(r)), and A£[q(r)] = *(g(t)) - 
$(q(0)), respectively [27[. They satisfy the energy con- 
servation AS = W — Q. Among these, the PDF for work 
P W {W) = (S(W[q(T)} - W)) satisfies the FT Q 



P W (W) 

p w {-w) 



(3) 



Recently, it was found that joint probabilities for ther- 
modynamic quantities also satisfy similar fluctuation re- 
lations [28| . Let {_Ai[<7(Y)]} be a set of functionals whose 
sum is equal to the work W[q(T)] = J^. A;[q(r)] and 
-4i[9( r )] = ~-4«[9( T )] where q(r) = q(t — t) is a time- 
reversed path. Then, it was found that the joint PDF 
P({Ai}) = (JIAMQ] ~ A i)) satisfies [H 



g(W) 

^({-4}) 



(4) 



Those relations reproduce Eq. <j3j> , and provide more de- 
tailed informations on nonequilibrium fluctuations [29I ] . 

We apply the formalism to the study of heat fluctua- 
tions. Consider the joint PDF P h . e (Q,AE) = (S(Q[q\ - 
Q)S(AS[q] - AE)). Since W = Q + AS, it follows from 
Eq. gj) that 



P h>e (Q,AE) = e Q+AE P h>e (-Q, -AE) 



(5) 



which will be referred to as a generalized FT (GFT). 
One may consider another joint PDF P w ^ e (W, AE) = 
(5{W[q] - W)5{AS[q] - AE)). They are related as 

P h , e {Q, AE) = P w , e {Q + AE, AE) . (6) 

So the GFT for P w . e takes a slightly different form as 

P w , e {W, AE) = e w P w , e {~W, -AE) . (7) 

The PDF P h (Q) = (S{Q[q] - Q)) for heat is reduced 
from Ph, e (Q, AE). Integrating both sides of Eq. (|5|) over 
(AE), we obtain a fluctuation relation for the heat: 

where 

*(Q) = J d(AE) e- AE P elh (AE\Q) . (9) 

Note that P e \ h (AE\Q) = P hie (Q, AE)/P h (Q) denotes a 
conditional probability for an energy change AE to a 
given value of heat dissipation Q. A reciprocity relation 
= ^(Q)" 1 was used in Eq. ©. The integral 
version is obtained from Eq. ([5} or Eq. ([5]) . It is given by 



e -S[q(r)]\ = / e -A£[q(r)] 



(10) 



The detailed FT for heat is modified by the factor 
The original FT requires that *(Q) = 1 for all Q. 
However one, in general, expects a correlation between 
Q and AE. Such a correlation leads to a Q-dependence 
in *f>(Q), hence invalidates the detailed FT for finite t. 

The fluctuation relations can be rewritten in terms of 
moment generating functions Gw(ty = (e~ AVV ), G^iv) = 
(e^Q), g^(X,K) = (e- xw - A£ ), and Q^k) = 
(e~ ??2_KA£ ). All of them are not independent, but are 
deduced from a single one, e.g., Gw,e- Equation (|6j) yields 
that S' hg (r], re) = Gw,e(Vi K ~ V)> and SwW = Gw,e{\ re = 
0) and '^(77) = G Lg (v^ = 0). The GFTs in Eqs. © 
and (0 arc equivalent to 



Gw,e(\ re) = Gw,e(l - A, -re) , 

Gt 6 (v, «) = Gt 6 (l - v, 1 - k) 



(11) 

(12) 



Setting re = in Eq. (fTTj) . one recovers the FT for work, 
GwW = GwO- — A). The fluctuation relation for heat can- 
not be written in a simple form with generating functions. 
Instead, the modification factor ^(Q) can be written as 

= QhAQ> K = V/GhAQ, re = o) (13) 

with Gh,e(Q, re) = / d(AE)e^^P Ke {Q, AE). 

In the t — > oo limit, the FT is formulated in terms 
of the large deviation function (LDF) For the heat 
distribution, it is defined as 



eh(q) 

Then, Eq. © yields that 



lim --hiP h (Q = qt) . 



where 



eh{-q) -eh(q) = q + i'(q) 



tfj(q) = lim ln*(Q = qt) . 



(14) 
(15) 
(16) 



We can further simplify it by introducing a LDF 

e h ,s(q,K)= lim -\ lnGh,e{Q = qt, re) , (17) 

t— >oo t 

which is obtained from the Legendre transformation 

eh,e(q, re) = max{e^ .(r), re) - qrj} (18) 

of a LDF e^ .(?7, re) = lim^oo - \ ]nGh,e(v> re) . Combin- 
ing these, we obtain that 



fpiq) = eh,e(q, re = 1) - e/,,g(g, re = 0) 



(19) 



This is a central quantity that determines whether the 
FT holds for heat in the f ^ oo limit. 

We apply the formalism to a d = 2 dimensional linear 
diffusion system where the force is given by f(q) = — f q 
with a force matrix 



1 e 
-e 1 



(20) 
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T\ q q T\ 



FIG. 1. Large deviation functions for the generation function 
in (a) and the PDF in (b) for heat at e — l/(2y / 3) are drawn 
with solid lines. Dashed lines are the plot of (fi(r?) — 1) in (a) 
and its Legendre transformation in (b). 



This model is a specific case of a general linear diffusion 
system studied in Ref. 271 ] . where one can find closed 
form solutions for various distribution functions. The 
purpose of this study is to confirm the fluctuation re- 
lations in Eqs. ([5]), Q, and © explicitly, and to un- 



derstand the effect of the correlation on the fluctuation 
relation. 

The force matrix is decomposed into the symmetric 
part F s = (F + F T )/2 = I and the anti-symmetric part 
F a = (F — F T )/2. Then, the conservative force, the en- 
ergy function, and the nonconservative force are given 
by f c (q) = ~q, #(g) = q T F s - q/2 = ±|qf, and 
f nc = — F a ■ q = e(—q2,qi) T , respectively. So the model 
describes a particle trapped in an isotropic harmonic po- 
tential $ and driven by a swirling force f nc . The param- 
eter e represents a strength of the driving force. 

The linear diffusion system was studied in Ref. [27| us- 
ing a path-integral formalism. We extend the formalism 
to obtain the joint probability distributions. The algebra 
is straightforward but rather lengthy. So we present the 
explicit expression for the moment generation function 
<?u>,e(A, re) without derivation. Details will be published 
elsewhere ■ It is given by 



Qw,s(\ k) = -^(A, re) 



(21) 



where 



F{x,y) 



cosh(f*(x)t) + (1 ~ 4 Sx" (a)2 aiah(n(g)t) 



and 



fi(x) = y/l - Ae 2 x(x - 1) . 



(22) 



(23) 



Note that f2(x) = f2(l — x). Hence, Gwe(\ K ) satisfies 
the GFT in Eq. (jTTj). 

The generating function for Ph(Q) is given by 



(24) 



FIG. 2. (a) Deviation from the FT for heat at e = 1/(2^3). 
(b) Validity domain for the LDF ~{r], re). 



It does not obey the FT (G h (v) ^ £7^(1 — 77)) for finite t. 
The corresponding LDF is given by 



e-iin- <V<V+) = tt(r)) - 1 



(25) 



We remark that the limit should be taken carefully. The 
function F(x,y) has a pole singularity at (f2(x) — l) 2 = 
Ay 2 in the t — > 00 limit. Hence, the LDF is well-defined 
only within the interval r\- < r\ < ?y+ where 77+ = 1 
and r}- = \— ^y/l + l/e 2 for e 2 > 1/3 and r/_ = (e 2 - 
l)/(e 2 + 1) for e 2 < 1/3. Equation (J25) is valid only 
within the interval, while e^(ry) = — 00 otherwise. 

The Legendre transformation eh(q) = max,{e,-(ij) — 
577} yields that 



eh(q) 



-r)+q + Sl(r}+) - 1 
— 7?_g + r2(?7_) - 1 



(l + E 2 )( 9 2 +4£^) _ g 

4e 2 2 



, <7 < <7+ 

, q>q- 

1 , otherwise 



(26) 



where q± = de~ h jdr\\^ . The linear branches indicate 
exponential tails in Ph (Q) . 

Figure (Ha) shows the LDF at e = l/(2\/3). The func- 
tion [0(77) — 1] is drawn with a dashed line, while e h (rj) 
is drawn with a solid line. The Legendre transformation 
eh (q) is plotted in Fig. QJb) with a solid line. The Leg- 
endre transformation of [Cl(rf) — 1] is also drawn with a 
dashed line. They deviate from each other at q±. 

In order to test the FT, we plot eh(—q) — &h{q) (solid 
line) in Fig. HJa). It does not coincide with the dashed 
straight line representing e^(— q) — eh{q) = q for large 
q, which shows that heat does not obey the FT. It is 
worthy to compare our result with that of Ref. [l9| • In 
both cases, the FT appears to be valid for small values of 
q, specifically within the domain \q\ < \q + \ in our study. 
The value of |e^(— q) — eh(q)\ saturates to a constant for 
large \q\ in Ref. [19(. It contrasts with the linear increase 
when |g| > \q_\ in our case. It suggests that the heat 
fluctuations do not exhibit a universal behavior 

Our theory predicts that ip{q) describes the deviation 
from the FT. We now evaluate tp(q) explicitly to confirm 
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the proposed relation in Eq. ([T5]) . Using Eq. (O and 
^M^' K ) = GwAv, K - V), one has Gj l & (v, = ^C??, « - 
7]). So the LDF is given by 

e M (r ?)K )=0(7 7 )-l. (27) 

It appears to be independent of k. However, due to the 
singularity of T in Eq. (f2"2"j) , the LDF is well defined only 
within the region 4(k — t?) 2 < (0(7?) — l) 2 , i.e., 

l + fifa) . 1 + 0(77) r9a s 

«7 j - K - ?/ 2 ' ^ 

Accordingly, the LDF ~(r], k) has a k dependence. The 
domain is drawn in Fig. [UJb). 

Now we need perform the Legendre transformation of 
Eq. (fTS)) at k — and 1. When n — 0, i) is restricted 
to the interval 77- < 77 < 77+, and eh,g{q, k — 0) becomes 
equal to eh(q) given in Eq. (|26|) . When k = 1, the validity 
region is shifted to 77^ < 77 < 77^ (see Fig. [2jb)). So, 
6h,e(q, K = 1) is given by the function in Eq. (f2l)| with 
7/± and g± being replaced with r)'± and q'± = de^/dr]^—^' , 
respectively. Notice the symmetry 0(77) = £1(1 — 77). It 
yields that 77^ = 1 — 77^ and q' ± = —q^- Inserting these 
into Eq. (|26p . one can find that eh,e{q, K = 1) = ej,(— g) — 
q. This completes the proof that ip(q) = &h[—q) — &h{q)~ 

q- 

We also test validity of the relation ([8]) at finite t. We 
have solved Eq. ([1]) with e = l/\/3 numerically 10 7 times 
up to t = 0.5 to measure various PDFs and ^(Q). In 
Fig.[3ja), P W (W) and e w P w (—W) are compared, which 
confirms the FT for work. In Fig. [^b) Ph{Q) displays 
a disagreement with e^Ph(-Q) but matches perfectly 
with Ph{—Q) l^{Q)- This is a numerical verification 
of the relation in Eq. ©. The joint PDF P w<e (W, AE) 
shown in Fig. [3tc) is symmetric under inversion AE — > 
—AE. So the energy may increase or decrease equally 
likely irrespective of the amount of work. It explains the 
reason why the heat distribution is wider than the work 
distribution as shown in Fig. [3] One can find an anti- 
correlation between Q and A£ in Fig. [3Jd). Due to the 
correlation, l'(Q) = (e~ AB ) Q ^ 1. 

In summary, we have derived the fluctuation relation 
for heat in Eq. flSJ using the GFT in Eq. ((5]) for the joint 
PDF. The heat distribution docs not obey the same type 
of the fluctuation relation as the work distribution does. 
The modification is given by the factor ^(Q) = (c~ ae )q 
that depends on the correlation between heat and energy 
change. The modified fluctuation relation for the heat 
has been tested analytically and numerically for a linear 
diffusion system. 

Our result shows that the FT for heat is not valid in 
general for finite t. The model studies in this work and 
in Ref. [3] show explicitly that the FT is violated even 
for the LDF in the t — > 00 limit. Nevertheless, it still 
remains as an open question whether there is a criterion 
for the FT in terms of the LDF. A sufficient condition is 




FIG. 3. (Color online) (a) Semi- log plots of P W (W) and 
e w P w (-W). (b) Semi-log plots of P h {Q), e Q P h (-Q), and 
e Q P h (-Q)/^(Q). Density plots of In P w , e (W, AE) in (c) and 
lnP h , e (Q,A) in (d). 



readily obtained from our result. Suppose that the en- 
ergy function is strictly bounded as Eq < $(q) < E\ with 
finite £0,1 HHIli- Then, e~^~ Eo) < *(Q) < e El ~ E \ 
hence ip(q) = and the FT holds. Hopefully, our formal- 
ism may yield a more strict condition for the FT. Future 
works are necessary in order to understand implication of 
the proposed fluctuation relation and to generalize it for 
systems with a time-dependent perturbation or systems 
in contact with many reservoirs. Experimental studies 
in small-sized systems [l6| are also necessary in order to 
characterize nonequilibrium fluctuations of heat. 
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